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Harmonic Sums, Polylogarithms, Special 
Numbers, and their Generalizations 

Jakob Ablinger and Johannes Bliimlein 

"en 

o 

(N 

U 

Abstract In these introductory lectures we discuss classes of presently known 
nested sums, associated iterated integrals, and special constants which hierarchi- 
cally appear in the evaluation of massless and massive Feynman diagrams at higher 
rS^ ' loops. These quantities are elements of stuffle and shuffle algebras implying alge- 

C^H, braic relations being widely independent of the special quantities considered. They 

r~| ' are supplemented by structural relations. The generalizations are given in terms of 

generalized harmonic sums, (generalized) cyclotomic sums, and sums containing 
in addition binomial and inverse-binomial weights. To all these quantities iterated 
integrals and special numbers are associated. We also discuss the analytic continu- 
ation of nested sums of different kind to complex values of the external summation 
bound N. 
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1 Introduction 



^^ ■ In the solution of physical problems very often new classes of special functions 

have been created during the last three centuries, cf . |QJ|5] . This applies especially 
also to the analytic calculation of Feynman-parameter integrals [6| for massless and 
massive two- and more-point functions, also containing local operator insertions 
and corresponding quantities, cf. [7, 8|. In case of zero mass-scale quantities the 
associated integrals map to special numbers, lately having been called periods |9|, 
Cd ■ see also [10|. In case of single-scale quantities, expressed as a ratio x G [0,1] to 

the defining mass scale, the integrals are Poincare iterated integrals 111 111121 or they 
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emerge as a Mellin-transform at N € N Ifl3l in terms of multiply nested sums. A 
systematic way to these structures has been described in I1141I15L Here an essential 
tool consists in representations by Mellin-Barnes [ 16] integrals. They are applicable 
also for integrals of multi-scale more-loop and multi-leg Feynman integrals [17|, 
which are, however, less explored at present. 

In the practical calculations dimensional regularization in D = 4 + e space-time 
dimensions ifTSl is used, which is essential to maintain conservation laws due to 
the Noether theorem and probability. It provides the singularities of the problem in 
terms of poles in £. However, the Feynman parameter integrals are not performed 
over rational integrands but hyperexponential ones. Thus one passes through higher 
transcendental functions J4]|5| from the beginning. The renormalization is carried 
out in the MS-scheme, chosen as the standard. In new calculations various ingredi- 
ents as anomalous dimensions and expansion coefficients of the j3 -functions needed 
in the renormalization can thus be used referring to results given in the literature. At 
higher orders the calculation of these quantities requests a major investment and is 
not easily repeated at present within other schemes in a short time. 

With growing complexity of the perturbative calculations in Quantum Field The- 
ories the functions emerging in integration and summation had to be systematized. 
While a series of massless 2-loop calculations, cf. [19], during the 1980ies and 90ies 
initially still could be performed referring to the classical polylogarithms 1 12 20-23 1 
and Nielsen-integrals [24|, the structure of the results became readily involved. In 
1998 a first general standard was introduced [25. 26| by the nested harmonic sums, 
and shortly after the harmonic polylogarithms lETl . Further extensions are given by 
the generalized harmonic sums, the so-called S-sums [28 , 29 1 and the (generalized) 
cyclotomic sums [30|, see Figure 1. Considering problems at even higher loops and 
a growing number of legs, also associated with more mass scales, one expects vari- 
ous new levels of generalization to emerge. In particular, also elliptic integrals will 
contribute [311. 




Fig. 1 Relations between the different extensions of harmonic sums. 
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These structures can be found systematically by applying symbolic summation, 
cf. [32), and integration formalisms, cf. 03311341 , which also allow to proof the rel- 
ative transcendence of the basis elements found and are therefore applied in the 
calculation of Feynman diagrams. 

In this survey we present an introduction to a series of well- studied structures 
which have been unraveled during the last years. The paper is organized as follows. 
In Section 2 a survey is given on polylogarithms, Nielsen integrals and harmonic 
polylogarithms. In Section 3 harmonic sums are discussed. Both harmonic polyloga- 
rithms and harmonic sums obey algebraic and structural relations on which a survey 
is given in Section 4. In Section 5 we discuss properties of the multiple zeta values 
which emerge as special constants in the context of harmonic sums and polyloga- 
rithms. The S-sums, associated iterated integrals, and special numbers are consid- 
ered in Section 6. The generalization of harmonic sums and S-sums to (generalized) 
cyclotomic sums, polylogarithms and numbers is given in Section 7. A further gen- 
eralization, which appears in massive multi-loop calculations, to nested binomial 
and inverse-binomial harmonic sums and polylogarithms is outlined in Section 8. 
Finally, we discuss in Section 9 the analytic continuation of the different kind of 
nested sums in the argument N to complex numbers, which is needed in various 
physical applications. Section 10 contains the conclusions. The various mathemat- 
ical relations between the different quantities being discussed in the present article 
are implemented in the package HarmonicSums . m 



2 Polylogarithms, Nielsen Integrals, Harmonic Polylogarithms 



Different particle propagators l/A<-(/7,,m,) can be linked using Feynman's integral 
representation 11361 

i _r(LUv k ) fft nLi^-; S (i-±x k ) ,v,eK.(i) 



While the momentum integrals over pi can be easily performed, the problem con- 
sists in integrating the Feynman parameters x^. In the simplest cases the associated 
integrand is a multi-rational function. In the first integrals one obtains multi-rational 
functions, but also logarithms l33l . The logarithms [37] have to be introduced as 
new functions being transcendental to the rational functions 

' L -ln(l-jc), etc. (2) 



/o 1 - z 
Iterating this integral by 
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dz\ /"« dz 2 , . ... 

— / i =Li2W (3) 

Zl JO 1— Z2 

one obtains the dilogarithm or Spence-function [20, 21], which may be extended to 
the classical polylogarithms [ 12 , 21-23] 

f x dz 

/ -Li n _ l (z)=U n (x),neN. (4) 

Jo z 

All these functions are transcendental to the former ones. For an early occurrence 
of the dilogarithm in Quantum Field Theory see | 



The above iterations are special cases in iterating differential forms in {dz/z, dz/ ( 1 
z)}. The general case is described as Nielsen integrals [24]. 

(•_iy»+p-l t \ dz 
S nj ,(x)= \ > , / ^ln- l (z)ln"(l- X z). (5) 

(n- \)\p\ Jo z 

Likewise, one might also consider the set {dz/z,dz/(l +z)}- Nielsen integrals obey 
the relation 

%-i,p(x) = -r%,p( x ) ■ ( 6 ) 

ax 

One may derive serial representations around x — 0, as e.g.: 

Li «M=£^> S 1 , 2 W = £^5 1 (fe-l), S 2 , 2 W=£^5 1 (fe-l), (7) 

see also 11391 . Here Si (ft) = L/Li(l/&) denotes the harmonic sum. The Nielsen in- 
tegrals obey various relations [12. 20-24]. A few examples are : 

Li 2 (l-x) = -Li 2 (x)-ln(x)ln(l-x) + C2 (8) 

Li 2 f-i) = -Li 2 (-x) - iln 2 (x) - C2 (9) 

Li 3 (l -*) = S l>2 (x) -ln(l -x)Li 2 (x) - iln(x)ln 2 (l -x) + £ 2 ln(l -x) + & 

(10) 
Li 4 I -—— J = ln(l -x)[Li 3 (x) - Si 2 (x)] + S 2 . 2 (x) -Li 4 (x) - Si, 3 (x) 



-V(l-x)Li 2 (x)-^ln 4 (l-x) (11) 

Li„(x 2 ) = 2"- 1 [U n (x)+U n (-x)} (12) 

Li 2 (z) = -£Li 2 (x),«eN\{0} (") 

S2,2(l -*) = -S2,2(x)+ln(x)Si 2 W- [Li 3 (x)-ln(x)Li 2 (x)- C 3 ]ln(l -x) 
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+ iln 2 (x)ln 2 (l-x) + ^. (14) 

Here £„ = Y,k=\ (1A*)> n — 2, n € N are values of Riemann's ^-function. 

Going to higher orders in perturbation theory it turns out that the Nielsen inte- 
grals are sufficient for massless and some massive two-loop problems, cf. M26II401 . 
and as well for the 3 -loop anomalous dimensions [41 1, allowing for some extended 
arguments as — x,x 2 . At a given level of complexity, however, one has to refer to a 
more general alphabet, namely 

21 = {cOo,Oh,CO-i} = {dz/z,dz/{l~z) 7 dz/{l +z)} . (15) 

The corresponding iterated integrals are called harmonic polylogarithms (HPLs) 
ll27ll . Possibly the first new integral is 

H_i,o,o,iW= / t~, — ■ (16) 

Jo z l+z 

Here we use a systematic notion defining the Poincare iterated integrals lHTl[T2l . 
unlike the case in (0. The weight w =\ HPLs are 

Ho(*)=ln(jc), Hi(jc) = -ln(l-jc), H_i(jc) =ln(l +x), (17) 

with the definition of Ho q(x) =ln"(x)/nl for all x indices equal to zero. The above 
functions have the following representation 

rx rx rx k 

Li„(x)=/ (Oq~ 1 (Oi, S p>tt (x)= co'qCo", U mw {x)= I Y\co m , , (18) 

JO JO JO 7_i 

where the corresponding products are non-commutative, m v , is of length k and x > 

Z\ > ... >Z m - 

Harmonic polylogarithms obey algebraic and structural relations, which will be 
discussed in Section|4] Numerical representations of HPLs were given in 114211431 . 



3 Harmonic Sums 



The harmonic sums are recursively defined by 

jv ( s \en(b)) k 
SbA N ) = E ,|L S a (k) , S«(N) = 1 , b,a t € Z\{0}. (19) 

k=\ k 
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In physics applications they appeared early in J44]|45]. Their systematic use dates 
back to Refs. [25, 26|. They can be represented as a Mellin transform 



S a (N) = M [/(*)] (N) = f dx ^- l /(*), N € N\{0}, (20) 

Jo 

where f(x) denotes a linear combination of HPLs. For example, 



24 4 

71n(2)C 3 , & 



(21) 



holds. Harmonic sums possess algebraic and structural relations, cf. Sect.|4] In the 
limit N — >• °° they define the multiple zeta values, cf. Sect. [5] They are originally de- 
fined at integer argument N. In physical applications they emerge in the context of 
the light-cone expansion |46|. The corresponding operator matrix elements are ana- 
lytically continued to complex values of Af either from the even or the odd integers, 
cf. Sect. H] 



4 Algebraic and Structural Relations 
4.1 Algebraic Relations 



Algebraic relations of harmonic polylogarithms and harmonic sums, respectively, 
are implied by their products and depend on their index structure only, i.e. they 
are a consequence of the associated shuffle or quasi-shuffle (sniffle) algebras [47|. 
These properties are widely independent of the specific realization of these algebras. 
To one of us (JB) it appeared as a striking surprise, when finding the determinant- 
formula for harmonic sums of equal argument [26 1 Eqs. (157,158) 

1 * 
Sa,...,a(A0=£LS a ,... )a (A0S A M a (A0, a f\ b = S ign(ab) (\a\ + \b\) (22) 

k ~ / 



also in Ramanujan's notebook H48II . but for integer sums, which clearly differ in 
value from the former ones. Related relations to again different quantities were given 
by Faa die Bruno [49|. 



Title Suppressed Due to Excessive Length 7 

Iterated integrals with the same argument x obey shuffle relations w.r.t. their prod- 
uct, 

H ai ,...,a k {x) •'Hb l ,...,bi{x) = Y n c u ...c k+l {x). (23) 

ceaUUb 

The shuffle-operation runs over all combinations of the sets a and b leaving the order 
of these sets unchanged. Likewise, the (generalized) harmonic sums obey quasi- 
shuffle or stuffle-relations, which are found recursively using [28 29 1 

S<n,...0 k (xi ,■■■ ,x k ;n) S h[ &/ (yi ,. . . ,yf,n) = 

n yi 

Y ^7 s «2v :«t( x 2, ■ ■ ■ ,**;») s &i b,(yu- ■ ■ ,yi,i) 

i=i l 

' y\ 



+ Y *7 s »i >-.o*( Jf i ' • • • >■**> s h b,(y2, ■■■,yv,i) 

v (xi-yi)' , ,v„ , .^ 

_ 2- -ai+fci ^...^{xz, ■ ■ ■ ,x k ,i) ^b 2 ,...,b,yy2, ■ ■ ■ ,yi,i) , 

jC(,y,eC,a/,t<GN\{0}. (24) 

The presence of trace terms in form of lower weight products in addition to the 
shuffled terms, cf. Il50l . leads to the name sniffle relations. In case the corresponding 
values exist, both (123124b can be applied to the multiple zeta values or other special 
numbers applying the integral and sum-representations at x — 1 and Af — >• °°, cf. [51 1. 
The basis elements applying the (quasi) shuffle relations in case of the harmonic 
sums and polylogarithms at a given weight w can be identified by the Lyndon words 
B52l53L Letgl-: {a,Z?,c,flf,...}beanorafere(ialphabetand2l*(2l) the set of words w 
given as concatenation products. Under the ordering of 21 a Lyndon word is smaller 
than any of its suffixes. For example, the set {a,a,a,b,b,b}, a < b is associated to 
the Lyndon words {aaabbb.aababb^aabbab}. Radford showed [54] that a shuffle 
algebra is freely generated by the Lyndon words. The number of Lyndon words can 
be counted using Witt formulae [55 1. Let 9Jt be a set of letters q in which the letter 
a k emerges n k times, and n = Ef =1 n k . The number of Lyndon words associated to 
this set is given by 

Unu...,n q )= 1 -Ym ( „, / J:/t / ^ . (25) 



// 



d\tik 



{ ni /d)\...{n q /d)\ 



Similarly one may count the basis elements occurring for all combinations at a given 
weight, if the alphabet has m letters : 

N A ( W ) = -Y^(^)'n d , (26) 

d\w 
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where pi denotes Mobius' function l56l . In case of the harmonic sums and polylog- 
arithms one has m = 3. The original number of harmonic polylogarithms is 3 W and 
in case of the harmonic sums 2 • 3 W ~ ' . Algebraic relations for the harmonic polylog- 
arithms and harmonic sums are implemented in the FORM-codes summer (25| and 
harmpol J27], HPL ll57l, and also HarmonicSums .m H29II35I . 



4.2 Structural Relations 



Structural relations of harmonic polylogarithms and harmonic sums are implied by 
operations on their arguments x and N, respectively. 



4.2.1 Harmonic Polylogarithms 



Harmonic polylogarithms satisfy argument-relations, as has been illustrated in (|8]- 
[Pfl i for some examples in case of the Nielsen integrals. Not all argument relations 
map inside the harmonic polylogarithms, however, cf. [27]. Some of them are valid 
only for the sub-alphabet {too, ©i }■ While the transformation x — > —x is general 

Ha(-x) = (-l) p H_ a (x), (27) 

with the last letter in a diffeent from and p the number of non-zero letters in a. 
The transformations 

x^l—x, X^rX 2 (28) 

apply to subsets only. Examples are : 

Hi,o,i(l-x) = -Ho(x)Ho.i(x)+2Ho.o.iM-C2Ho(x)-2C3 (29) 

Hi A o,i (* 2 ) =4[H 1| o,o,i(*)-H 1 ,o,o,-i(*)-H_i A o,i(jf)+H_i i o f0 _i(*)].(30) 

One may transform arguments by x — >• \/y + '£, 

H i,0,i(~) =H (x)[Ho i i(x)+^H 1 (x)-4C2 + 7r 2 ]-2[H OA1 (x)-Ho,i,i(x) 

+ C 3 ] + [-Hi(x) -»7t]Ho,i(jO +2£tHi(x) - U\l{x) 

o 

+~i[7C + iHi(x)}Kl(x) . (31) 



Title Suppressed Due to Excessive Length 
An important general transformation is 



x— > 



\-t 

T+7 



(32) 



which acts on the HPLs but not on the subset of Nielsen-integrals. An example is : 

Hl, ~ 1,0 (tTjc) = J H -iW+ H -i-i.iW- H o,-i,-iW-Ho,-i,iW + ^ 



1 



&p_i(je)-Ho(.x)]-2 



C3 ln(2)C 2 



-21n(2)C 2 



(33) 



In most of these relations also HPLs at arguments = 1 contribute, cf. Sect. [5] Struc- 
tural relations of HPLs are implemented in the packages harmpol ll27l . HPL (53, 

and HarmonicSums .m I 



4.2.2 HarmonicSums 



Harmonic sums obey the duplication relation 

S iu ...,i n {N) = 2"+- + '"-"£5 ± , b ..., ±I „(2^V), i k e N\{0} . (34) 

± 

This allows to define harmonic sums at half-integer, i.e. rational, values. Ultimately, 
one would like to derive expressions for JVeC, cf. Sect. [9] Another extension is 
to N G M 1 14 26 1 . The representation of harmonic sums through Mellin-transforms 
d20b implies analyticity for a finite range around a given value of N. The Mellin- 
transform of a harmonic polylogarithm can thus be differentiated for N 



d f 1 f l 

— / dx^-^Jx) = / dxx N - 1 H (x)H a (x) . 
dN Jo Jo 



(35) 



In turn, the shuffling relation ( 123) allows to represent the r.h.s. in d35T > as the Mellin- 
transform of other HPLs. It turns out that differentiation of harmonic sums for JV is 
closed under additional association of the multiple zeta values 1 14 1. The number of 
basis elements by applying the duplication relation (H), resp. its combination with 
the algebraic relations is 



N H (w) =2-3 w ~ 1 -2 w -\ %W = -£|i- 



w • 



2 2 _ yl 



(36) 



Differentiation in combination with the other relation yields 
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N D (w)=4-3 w -\N DH (w)=4-3 w - 2 -2 w - 2 ,N ADH {w) =N ah {w) -N ah (w- 1). 

(37) 



Let us close with a remark on observables or related quantities in physics which 
are calculated to a certain loop level and can be thoroughly expressed in terms of 
harmonic sums. As a detailed investigation of massless and single mass 2-loop quan- 
tities showed |j40ll seven basic functions of up to weight w = 4, cf. Ifl4l . are suffi- 
cient to express all quantities. The 3-loop anomalous dimensions BD contributing 
to the 1/e poles of the corresponding matrix elements require fifteen functions of 
up to weight w = 5 and further twenty basic functions are needed to also express 
the massless Wilson coefficients [7] in deep-inelastic scattering [59], cf. Ref. [60|. 
Despite of the complexity of these calculations finally a rather compact structure is 
obtained for the representation of the results. Structural relations of harmonic sums 
are implemented in the package HarmonicSums . m | 



5 Multiple Zeta Values 



The multiple zeta values (MZVs) [61 62 Q| are obtained by the limit N — > °° of the 
harmonic sums 



lim S B (JV) = da (38) 



N 

and may also be represented in terms of linear combinations of harmonic polylog- 
arithms Hb(l) over the alphabet {(Bo,a)i,0)_i}o In the former case one usually 
includes the divergent harmonic sums since all divergent contributions are uniquely 
represented in terms of polynomials in <7i (°°) = <j due to the algebraic relations. 
Likewise, not all harmonic poly logarithms can be calculated at x — 1, requiring their 
re-definition in terms of distributions. Some examples for MZVs, which already ap- 
pear in case of Nielsen integrals, are : 

Li„(l) = ;„ (39) 

1 
2" 
Si, P (l) - Cp+1 (41) 



U,(-l)=-(l-— )C» (40) 



1 For a detailed account on the literature on MZVs see [63. 64 1 and the surveys Ref. |65 1. 

2 The numbers associated with this alphabet are sometimes also called Euler-Zagier values and 



those of the sub-alphabet {cOq, a>\ } multiple zeta values. 
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Si, 2 (-1) = ~6 (42) 

Li 2 (£)=i[C2-ln 2 (2)] (43) 

Li 3 Q) = 1 -^- l -^H2)+ l -\n\2) (44) 

S 2 , 2 (l) = ^C 2 2 (45) 

S 2 . 2 (-l) = -^ 2 +2Li 4 Q)+^3ln(2)-ic 2 ln 2 (2) + ^ln 4 (2) (46) 

Si. 3 (-1) = -^C 2 2 + Li 4 Q) +^3ln(2)-ic 2 ln 2 (2) + ^ln 4 (2) . (47) 

In case of physics applications, MZVs played a role in loop calculations rather early, 
cf . 11661 . Since for Li„, ( 1 /2 ) for m — 2 , 3 these numbers are not elementary, (1461 l47t 
seem to fail to provide a corresponding relation for m = 4. Similarly, for larger 
values of m also no reduction has been observed. 

A central question concerns the representation of harmonic sums in terms of 
polynomial bases. This has been analyzed systematically in [63 67 1. [£| For MZVs 
over {0, 1} a proof on the maximum of basis elements at fixed weight w has been 
given in Refs. [68]. At the lowest weights the shuffle and stuffle relations imply all 
relations for the MZVs. Starting with weight w — 8 one also needs the duplication 
relation (l34l , and from weight w — 12 also the generalized duplication relations 
Sect. 4.1 in [63 1 . The latter are closely related to the conformal transformation rela- 
tions of the HPLs at x = 1, see (l32l . Let us give one example for the combined use 
of the shuffle and stuffle relation for illustration, [51 Q : 

shuffle : £ 2 .iC 2 = 6£ 3 , u +2£ 2 , 2 ,i + £2,1,2 
stuffle : £2.1 C 2 = 2 C 2 , 2 ,i + Cu + C 2 .3 + £2,1,2 

=► C3.1.1 = ^[C4.1 + C 2 ,3-C 2 . 2 .l] • (48) 

6 

Finally one derives a basis for the MZVs using the above relations. Up to weight 
w = 7 reads, cf. [25|, 

|((T 1 H,ln(2));C 2 ;C3;Li4( 1 5);(C5,Li 5 (i));(Li 6 (I) I (T_ 5 ,_ 1 ); 
($7 > Ii 7 (£),CT_ 5iM ,a_5 ) i,i);l (49) 



For some aspects of the earlier development including results by the Leuven-group, Zagier, 
Broadhurst, Vermaseren and the Lille-group, see 11631 . 
4 Here the £„ -values are defined £,, ,.„ m = £^ >„,>.„>„,„ UT=i n k "' 
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In ||63l bases were calculated up to w = 12 for the alphabet {0, 1,-1} and to w = 22 
for the alphabet {0, 1} in explicit form resp. for w — 24 restricting to basis elements 
only. In the latter case the conjecture by Zagier [62 1 that the shuffle, stuffle and 
duplication relations were the only ones was confirmed up to the weights quoted. For 
these cases counting relations were conjectured in Refs. [69 70|. One may represent 
the basis for the MZVs over the alphabets {0,1} resp. {0,1,-1} by polynomial 
bases or count just the factors appearing in these polynomials of special numbers 
occurring newly in the corresponding weight, which is called Lyndon-basis [63 j. 
In the first case the basis for the MZVs[{0, 1}] is conjectured to be counted by the 
Padovan numbers /\ [ 7 1 1 generated by 

i 1 2^ 3 = 1X4 P1=P2=P3 = 1- (50) 

1 X X k=Q 

In case of the Lyndon basis the Perrin numbers P^ appear ll72l 

o 2 °° 

,1* 3 = L**A, Pi=0,P 2 =2,P 3 ,= 3. (51) 

1 x x k=0 

Both the above sequences obey the Fibonacci-recurrence IT731 

P d =P d -2+Pd-3, d>3. (52) 

The length of the Lyndon basis at weight w is given by 

d\w 

Hoffman ll74l conjectured that all MZVs over the alphabet {0,1} can be represented 
over a basis of MZVs carrying 2 and 3 as indices only. This has been confirmed up 
to w = 24. An explicit proof has been given in ll75l . 

The polynomial basis of the MZVs[{0, 1, — 1}] is conjectured to be counted by 
the Fibonacci numbers ll76l 



f 1 

fd ~vs 



/-\ d / r-\ d ' 

1+V5\ / 1-V5\ 



(54) 



which obey 



1 x x k=Q 
For the corresponding Lyndon basis the counting relation 



= £**/*. (55) 
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l(w) = -l t ll(~)L d (56) 

d\w 

is conjectured 16911701 , where L c / are the Lucas numbers 11771 , 

2 V £**!*, (58) 



* X X k=0 



L d =L d _i+L d _ 2 ,d>4, L 1 = 1,L 2 =3,L 3 =4. (59) 



There is a series of Theorems proven on the MZVs, see also 11631 , which can be 
verified using the data base ll63l . The duality theorem 11621 in case of the alphabet 

{(Oq,(Oi} states 

H a (l)=H a1 (l), a^a^i. (60) 

In case of the alphabet {<Bo,a)i,0)_i} it is implied by the transformation (|32t , see 
ll63l. Another relation is the sum theorem, Ref. 



E &....A = C« ■ (61) 

i\+...+ik=n,i\>\ 

The sum-theorem was conjectured in ll80l . cf. JST]. For its derivation using the Euler 
connection formula for polylogarithms, cf. 



Further identities are given by the derivation theorem, Q80l[83J. Let / = 

(i i , . . . , 4) any sequence of positive integers with i\ > 1 . Its derivation D(I) is given 
by 

D (I) = (h + 1, h, ■ ■ ■ , h) + (h, h + 1) • • • , h) + ■ ■ ■ (h ,h, ■ ■ ■ , h + 1) 

Cd(/) = C(i l+ u 2 ,...,i k ) + ■■■ + C(i,i 2 i k +i) ■ (62) 

The derivation theorem states 

Cd{i) = Ct(d(t(/))) i (63) 

where T denotes the duality-operation, cf. d60l l. An index-word w is called admissi- 
ble, if its first letter is not 1 . The words form the set $j°. \w\ = W is the weight and 
d(w) the depth of w. For the MZVs the words w are build in terms of concatenation 
products Xq 1 x\Xq xq...x 1 q x\. The height of a word, ht(vv), counts the number 
of (non-commutative) factors Xqx\ of w. The operator D and its dual D act as fol- 
lows El, _ _ 

Dxq = 0, Dx\ — xqX\ , Dxq — xqx\ , Dx\ = . 
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Define an anti-symmetric derivation 

<9„.*o=-*o(-*o+-*i)"~ 1 *i • 

A generalization of the derivation theorem was given in [83. 85 1. The identity 

C(<9„w)=0 (64) 

holds for any n > 1 and any word w £ $}°. Further theorems are the Le— Murakami 
theorem, [86], the Ohno theorem, (83, which generalizes the sum- and duality 
theorem, the Ohno-Zagier theorem, [88], which covers the Le-Murakami theo- 
rem and the sum theorem, which generalizes a theorem by Hoffman [80 81 1, and 
the cyclic sum theorem, |89|. 



There are also relations for MZVs at repeated arguments, cf. 
which examples are : 



C({2}„) 

(2,{1}„) 

C({3,1}„) 



2(27t 



\2n 



i 



(2» + l)!2 
C(2,{l}„) = C(» + 2) 



4^ {i in) ~ (4n + 2)\ 



C({io}„) = 



In 



10(2tt) 
(10n + 5)! 



1+VB" 



10n+5 



1-V5 



10h+5' 



, on 

(65) 
(66) 
(67) 

(68) 



Finally, we mention a main conjecture for the MZVs over {0, 1 }. Consider tuples 
k= (k\,...,k r ) £ W,k x > 1. One defines 



■-■£\V 



{0} 

E Q-C(k)c 



|k|=vv 



(69) 



If further 



<.Go 



"Ca 



£ 2? w C E (Goncharov) 

:= 0JT w (Cartier) 

iv=0 



(70) 
(71) 
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the conjecture states 

(a) 2f Q ° = 2? c& . There are no relations over <Q> between the MZVs of different 
weight w. 

(b) dim J^, = d w , with do = 1 ,d\ = 0,^2 = 1 > d w = d w -2 + d w -3- 

(c) All algebraic relations between MZVs are given by the extended double-shuffle 
relations [91], cf. also ||92l . If this conjecture turns out to be true all MZVs are 
irrational numbers. 

Let us also mention a few interesting relations for Li2 (z) for special arguments 
found by Ramanujan |93|, which are of use e.g. in massive calculations at 3- 
loops [94|. These numbers relate to constants beyond the MZVs, which occur for 
generalized sums and their extension allowing for binomial and inverse binomial 
weights : 

Li2 (~£) + o" Li2 G) =-^ +ln ( 2 )M 3 )-^ ln2 ( 2 )-5 ln2 ( 3 ) C73) 
Li2 (l) + l U2 (l) =^+ 21n ( 2 ) ln (3)-21n 2 (2)-^ln 2 (3) (74) 

■*(-!)>(l)-5 + ^ 



L H-8j + 3 L H9j = -2 ln UJ- (?6) 

For further specific numbers, which occur in the context of Quantum Field Theory 
calculations see also Sections |6*31 171 andj9] 



6 Generalized Harmonic Sums and Polylogarithms 



6.1 Generalized Harmonic Sums 



Generalized harmonic sums, also called S-sums, are defined by 

N x h 

Sai a k {xi,---,Xk',N) = 2_i —Sa 2 ,...,a k ( x 2,---,Xk',il), 

'1=1 'l 
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5 = 1 , Xi e K\{0}, at £ N\{0} (77) 

and form a quasi-shuffie and a Hopf algebra [96] under the multiplication (124-b J28ll . 
The multiplication relation in general leads outside the weight sets {«,■},{&,•}. The 
S-sums cover (together with the limit N — > °o) the classical polylogarithms, the 
Nielsen functions, the harmonic polylogarithms, the multiple polylogarithms [97 1, 
the two-dimensional HPLs (98), and the MZVs El. In Ref. EH) four algorithms 
were presented allowing to perform the £-expansion of classes of sums in terms of 
S-sums, which were coded in two packages 09911 1001 . In this way the £-expansion 
can be performed using convergent serial representations for the generalized hy- 
pergeometric functions pFg, The Appell-functions Fi 2 > an d tne Kampe de Feriet 
function 11011 . 

They can be represented in terms of a Mellin transformation over x E [0,*i.. JtjJ 
||29l . E.g. the single sums are given by 



f b dx,„ [ x i dx 2 f X2 , jcf - 1 
Generalized harmonic sums obey the duplication relation 



JO X m Jo X2 Jo xi — l 



(78) 



£s flm ,..., a ,(±fc m ,-,±fc 1 ;2 jy) = 2l J ai _ m s am ,... fil (bl,...,b 2 l ;N) , (79) 

where the sum on the left hand side is over the 2'" possible combinations concerning 
± and at G N, bi € K\{0} and n £ N. They also obey differential relations w.r.t. 
N, supplementing their set with the generalized harmonic sums at infinity, resp. of 
the generalized harmonic polylogarithms at X = 1 . The mapping will usually also 
require objects with different weights x,-. Examples are 



3-S 2 (2;n) = -S 3 (2;n) +H (2)S 2 (2;n) +H ,o -i(l) +2Ho, ,i(l) +H ,i _i(l) , 
an 



^S 3 (i;A0 = 12 



-S 3> i(i, i;JV) - I As 2>1 ( i , i ;JV) - I Hl , (I) S 2 (I;A0 



U. m„ /,s I 



+Hq( 2 ) S 2 ,i(3, 2 ;iV) - -Hi ( 3 ) H 0i i,o( 2 ) + t^Ho.o.i.o 



2"^' ^'"^y 12 - 

i ,. i u. /, ... i. 



YWi)"^ 53 ^"^^ 



(80) 



The counting relations for the basis elements are 

N D (w)=N s {w)-N s {w-l), N A n{w)=N A (w)-N A (w-l), (81) 
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where Ns = (n — 1) -rC denotes the number of sums, given n letters in the alphabet, 
and Na the basis elements after applying the algebraic equations. Explicit bases for 
a series of alphabets have been calculated in |29|. 



6.2 Generalized Harmonic Polylogarithms 



Generalized harmonic polylogarithms are defined as the Poincare-iterated integrals 

MM 

rx m j 

h *w = / n i— i — ■ J ^ > °j g c ' z j ^ z /+i ■ ^ 82 > 

Jo + = i \a.j\ -sign(a j)Zj 



For aj e R, aj < 1 ( |82] | is defined as Cauchy principal value only. Already A. Jon- 
quiere lTT2l has studied integrals of this type. Sometimes they are also called Chen- 
iterated integrals, cf. [ 1 1 1, or Goncharov polylogarithms [97 ]. 

The Mellin transforms of generalized harmonic polylogarithms map onto gener- 
alized harmonic sums [29|. Furthermore, the generalized harmonic polylogarithms 
obey various argument relations similar to the case if the HPLs, cf. Sect. 14.2.11 as 

(83) 
(84) 

(85) 

(86) 

(87) 



6.3 Relations between S-Sums at Infinity 



S-sums at infinity exhibit a more divergent behaviour than harmonic sums if a\ > 1 . 
The degree of divergence is then at least °< a(, cf. Sect. [9] In the following we will 
discuss only convergent S-sums at infinity. They obey stuffle and shuffle relations, 
the duplication relation Af — > 2N, and the duality relations for the generalized poly- 
logarithms 



x + b 


— ► X 


b — x — > x 


1-x 


-> X 


l+x 


kx 


— > X 


1 


-* X 
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l-x->x (88) 

T+f -> x ( 89 ) 

gf -*•*, c,deR,d^l . (90) 



Eq. ( T88l > implies 



H aij ..., ajt (l)=H 1 _ fll> ..., 1 _ at (l), a*^0. (91) 

Examples for ([89]|90]l are : 

5i(i;oo)=-5_ 1 H = ln(2) (92) 

Si(f,oo)=-5_iH+5i(-i;oo) . (93) 

In various cases S-sums at infinity reduce to MZVs, cf. also [ 102|, 

S 1 ,i,i(i ) 2,l;°°)=~C2ln(2) + ~C3 (94) 

5 2 ,i(^l;-)=-^ln(2)C 2 + C3 (95) 

S M (^;oo)=Ii m (i) . (96) 

Otherwise, new basis elements occur which have both representations in infinite 
sums and iterated integrals. Using the above relations bases for different sets of 
S-sums at infinity were calculated in [ 1 02 1 . 



7 Cyclotomic Harmonic Sums and Polylogarithms and their 
Generalization 



The alphabet of the harmonic polylogarithms (Q3) contains two differential forms 
with denominators, which form the first two cyclotomic polynomials: (1 — x) and 
(1 +x). It turns out that quantum field theoretic calculations are also related to cy- 
clotomic harmonic polylogarithms and sums 01031 . Cyclotomic polynomials are de- 
fined by 

*n to = r?~\ ; > d , n G N+ (97) 

d\n.d<n 



and the generating alphabet reads 
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21= — ^U 



dx 1 f x dx 



k€N+,0< I <(p(k)\, (98) 



x J [ #*(*) 

where 4>fc(jc) denotes the fcth cyclotomic polynomial IU04L and (p(k) denotes Eu- 
ler's totient function 01051 . The Poincare iterated integrals over the alphabet d98l 
are called cyclotomic harmonic polylogarithms, cf. ll30l . Due to the regularity of 
l/<P n (x) for x G [0,1], except for <P\ (x), no more singularities appear beyond those 
known in the case of the usual harmonic polylogarithms (or Nielsen integrals). Cy- 
clotomic harmonic polylogarithms obey shuffle relations, cf. Sect. [4] 

The cyclotomic harmonic sums [30] are related to the cyclotomic harmonic poly- 
logarithms via a Mellin transform (l20V The generalized cyclotomic harmonic sums 
are given by 

S{ ai ,bl,ci} {a h b h ci}( s U---,sr-, N ) = 

N s k 

? ^ aikl l bl )c l S {"2.b2,c 2 };-;{''i^ h c l }^2,-,sr,ki) : S @ = l : (99) 

where ai,Ci G N+,£>; G N, si G R\{0},a,- > bf, the weight of this sum is defined 

by c\ H +c; and {a,-,fo/,c,} denote lists, not sets. If s, = ±1 these are the usual 

cyclotomic harmonic sums. The simplest cyclotomic sums are the single sums 

i.e. harmonic sums with cyclic gaps in the summation. The cyclotomic harmonic 
sums obey quasi-shuffle relations (A). 

Beyond this the cyclotomic harmonic sums obey structural relations implied by 
differentiation for the upper summation bound N, (D), which require to also con- 
sider their values at N — > °°. There are, furthermore, multiple argument relations, 
cf. 11301 , decomposing S aj j, h ci(k'N), called synchronization (M), and two duplica- 
tion relations (H\,H2). Let us consider the cyclotomic harmonic sums implied by 
the letters 

1 (-1)* 1 (-If 

— - — — — - — - — (101) 

#1' k'2 ' (2/t+l)'3' (2fc+l)'4 ' ^ ' 

The length of the basis can be calculated by 

N s (w) =4-5 n '- 1 (102) 

d\w 

Nd{w)=N s (w)-N s (w-1) (104) 
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»a,dW - N A {w)-N A (w- 1) (105) 

iW*M-iEMG)C#-3.20--l I E„(^W-3-* 



w rr \a / w- * , , 

aw rf w— 1 



i^ 



(106) 



where Ns(w) denotes the number of all sums. One may calculate the asymptotic 
representation of the cyclotomic harmonic sums analytically. Here also the values 
of cyclotomic sums at N —¥ °° occur. The singularities of the cyclotomic harmonic 
sums with Sk = ±1 are situated at the non-positive integers. 

The cyclotomic sums for N — > °° are denoted by 0{ ai bl C1 j {a h bi.ci}( s U---7 s i)- 
For V |^| < 1 divergent series do only occur if a\ = b\ and c\ — 1, where the de- 
gree of divergence is given by Ob as in the case of the harmonic sums and can be 
represented algebraically. They are related to the values of the cyclotomic harmonic 
polylogarithms at x = 1 . At w = 1 the regularized sums may be given in terms of 
\\f{k/l) and for higher weights in terms of \j/^"''(k/l),m > 1. If / is an integer for 
which the Z-polygon is constructible one obtains representations in terms of alge- 
braic numbers and logarithms of algebraic numbers, as well as K [30). In this way, 
£2 being a basis element in case of the MZVs, looses its role. At depth w — 2 Cata- 
lan's constant [106] with 

02,1,-2 = ~1+C (107) 

contributes. At higher depth new numbers emerge, which partly can be given inte- 
gral representations involving polylogarithms and roots of the integration variable 
x. The cyclotomic sums at infinity, as real representations, are closely related to 
the infinite generalized harmonic sums at weights Sk which are roots of unity, cf. 
also 111071 . In POl basis representations were worked out for w = 1,2 for the /th 
roots, / G [ 1 , 20] , cf . also [ 1 08 1 . Counting relations for bases of the cyclotomic sums 
at infinity have also been derived in Ref. [30|. 



8 Nested Binomial and Inverse-binomial Harmonic Sums and 
Associated Polylogarithms 



In massive calculations further extensions to the nested sums and iterated integrals 
being discussed in the previous sections occur. Here summation terms of the kind 
5 a .b,c( x ^)> where d = {di, ...,d m }, or their linear combinations are modulated by 

S a ,b.c(x;£) -> f jS a , b , c (x;fc) 
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S a ,b,c(x;£) -► -j^S^frk) , (108) 

k 

building iterated sums [94|. Sums of this kind occur in case of V-type 3-loop graphs 
for massive operator matrix elements. Simpler sums are obtained in case of 3-loop 
graphs with two fermionic lines of equal mass. Single sums of this kind have been 
considered earlier, see e.g. [ 109 ]. One may envisage generalizations of ( 1108b in 
choosing for the binomial a general hypergeometric term, i.e. a function, the ratio of 
which by all shifts of arguments being rational. The association of the correspond- 
ing iterated integrals in the foregoing cases has been found easily. Here the situation 
is more difficult and the functions representing these iterated sums are found in es- 
tablishing differential equations 111 01 . It is found in the cases occurring in Ref. 1194) 
that the corresponding differential equations finally factorize and one obtains iter- 
ated integrals over alphabets which also contain root-valued letters 

1 1 1 

a,b,c£Q (109) 



sj{x + a)(x + by ^(x + a){x + b)x + <. 

beyond those occurring in case of the generalized (cyclotomic) polylogarithms. A 
few examples of this type have been considered in 11111 . The relative transcendence 
of the nested sums and iterated integrals has been proven. The V-type 3-loop graphs 
require alphabets of about 30 root-valued letters. The corresponding nested sums do 
partly diverge « a N ,a eM,a>2. A typical example for a nested binomial sums is 
given by : 

E(T)(-2)'t-7^Mi.-i;J) (no) 



IV'V ' ^1 :N 



dx — r\/^[ H -i7,-i,oM- 2H ; i8 ,-i.oM] 

o x— 1 V 8 +x "' ' is. . j 



« / V ^ \ /^[Hf 2 W-2H? 3 (*)]+C3 /V 8V,V ' 
Jo x+1 \ 8+x Jo 



2 Jo x+1 VS+^ Jo x+i V 1 "*' 

with cj, = LJ=o^i.2 (ji - 1>7") U , -) 2 /j/(2j) ] -/ n one °f me specific constants emerg- 
ing in case of these sums. Here the iterated integrals H* extend to x = 1 as firm 
bound, contrary to the cases discussed before where x = is chosen. Here the new 
letters w^ are 

1 1 

W12 = —F= n VVB = 



y/x(8- x y ' (2-x)^x{%- x y 

w\-i = — , w\% = . (HI) 

y/x{%+x) ' (2+x)^/x(f+x) 
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The representations over the letters (1109) are needed to eliminate the power growth 
°c a N of these sums and can be used to derive the asymptotic representation at large 
values of N. While the terms °< 8 W and <* 4 N cancel, it may occur that individual 
scalar diagrams exhibit contributions °< 2 N , cf. (94). This behaviour is expected to 
cancel in the complete physics result. 



9 Analytic Continuation of Harmonic Sums 



The loop-corrections to various physical quantities take a particular simple form in 
Mellin-space being expressed in terms of harmonic sums and their generalizations. 
Moreover, in this representation the renormalization group equations can be solved 
analytically, cf. 159111 121 . For a wide variety of non-perturbative parton distributions 
Mellin-space representations can be given as well, see e.g. 11131 . 



Im(N) 




Re(N) 



Fig. 2 Path of the contour integral ( 11121 . 

Thus one obtains complete representations for observables in A^-space. In case of the 
perturbative part, the singularities are situated at the integers N < No, with usually 
Nq = I, see e.g. 01141 . The harmonic sums possess a unique polynomial represen- 
tation in terms of the sum Si (N) and harmonic sums which can be represented as 
Mellin transforms having a representation by factorial series 0115111161 . They are 
transformed to x-space by a single precise numerical contour integral around the 
singularities of the problem to compare with the data measured in experiment. The 
analytic continuation of the perturbative evolution kernels and Wilson coefficients 
from even or odd integers to complex values of N is unique [ 1 17 1£[ To perform 
this integral a representation of the harmonic sums for N G C is required. Accu- 
rate numeric representations up to w = 5 have been given in 11181 . see also [ 1 19 1. 



For a detailed proof also in case of generalized harmonic sums see 1291 
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Arbitrary precise representations can be obtained using the analytic expressions for 
the asymptotic representation [14 60] together with the recursion relations given in 
Sect. [3] These are given up to w — 8 in [58 1. A path to perform the inverse Mellin 
transform 

fix) = - I dz\m IV* x~ C M\f](N = CYl , C = c + ze^ (112) 

7t JO 

is shown in Figure 2. The asymptotic representations can also be obtained in ana- 
lytic form for the S-sums [291, cyclotomic (S)-sums [301, as well as f° r the nested 
binomial cyclotomic S-sums [94 110|. 

In case the expressions in TV-space result from Mellin transforms of functions 
f(x) ~ x a ,a e]0, 1[ the singularities are shifted by a. This is usually the case for 
the non-perturbative parton distribution functions, but also in case of some root- 
valued harmonic polylogarithms considered in Sect. [8] 

The inverse Mellin-transform cannot be performed in the above way for inte- 
grands which do not vanish sufficiently fast enough as \N\ — > °°. Contributions of this 
kind are those leading to distribution-valued terms in x-space as to 5(1 — x), [ln*(l — 
x)/(l — x)]+ the Dirac 8 -distribution and the +-distribution defined by 



/ dx [f(x)} + g(x)= [ dx [g(x)-g(l)]f(x) 
Jo Jo 



(113) 



Also for terms which grow like a N ,a G M.,a > 1 in TV-space, the Mellin transform 
cannot be performed numerically in general. They are not supposed to emerge in 
physical observables. The physical quantities in hadronic scattering contain the par- 
ton distribution functions, which, however, damp according contributions occurring 
in the evolution kernels sufficiently. On the other hand, the inverse Mellin transform 
can always be performed analytically changing form nested sum-representations in 
TV-space to iterated integral representations in x-space as has been outlined before. 



10 Conclusions 



Feynman integrals in Quantum Field Theories generate a hierarchic series of spe- 
cial functions, which allow their unique representation. They emerge in terms of 
special nested sums, iterated integrals and numbers. Their variety gradually extends 
enlarging the number of loops and legs, as well as the associated mass scales. The 
systematic exploration of these structures has been started about 15 years ago and 
several levels of complexity have been unraveled since. The relations of the various 
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integral representation (inv. Mellin transform) 




power series expansion 

Fig. 3 Connection between harmonic sums (H-Sums), S-sums (S-Sums) and cyclotomic harmonic 
sums (C-Sums), their values at infinity and harmonic polylogarithms (H-Logs), generalized har- 
monic polylogarithms (G-Logs) and cyclotomic harmonic polylogarithms (C-Logs) and their val- 
ues at special constants. 



associated sums and integrals are schematically illustrated in Fig. 3 and are widely 
explored. Many relations are implied by the shuffle resp. stuffle algebras, others are 
structural relations. The number of relations grows with the number of admissible 
operations. 

A large amount of transformations and relations between the different quantities 
being discussed in this article are encoded in the package HarmonicSums M29II35I 
for public use. For newly emerging structures the algebraic relations are easily gen- 
eralized but they will usually apply structural relations of a new kind. With the 
present programme revealing their strict (atomic) structure, they are fully explored 
analytically and Feynman's original approach to completely organize the calculation 
of observables in Quantum Field Theory is currently extended to massive calcula- 
tions at the 3-loop level in Quantum Electrodynamics and Quantum Chromodynam- 
ics at the perturbative side. For these quantities efficient numerical representations 
have to be derived. Working in Mellin space the treatment may even remain com- 
pletely analytic, in a very elegant way, up to a single final numerical contour integral 
around the singularities of the problem, cf. Sect. 



Despite of the achievements being obtained many more physical classes still 
await their systematic exploration in the future. It is clear, however, that the various 
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concrete structures are realized as combinations of words over certain alphabets, 
which may be called the genetic code of the microcosm [1201. 
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